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1 Introduction and summary 



Effective field theory allows one to analyze the chiral structure of Quantum Chromodynamics 
in the low energy domain, which is not accessible to a perturbative expansion in the strong 
coupling constant. The spontaneous violation of the chiral symmetry of QCD entails the 
existence of Goldstone bosons, the pions (we consider here the two flavor case of the light 
up and down quarks). The interactions of the Goldstone bosons among themselves and with 
matter vanish as the momentum transfer goes to zero. This is a consequence of Goldstone's 
theorem. Consequently, such interactions can be analyzed in a perturbative expansion, where 
all momenta and masses are small compared to the typical scale of hadronic interactions, say 
the mass of the rho meson. This is the so-called chiral expansion. A systematic investigation of 
processes involving pions allows therefore to understand in a precise and quantitative manner 
how the symmetry violation takes place and also to pin down the ratios of the light quark 
masses. One such process is elastic pion-nucleon scattering. In ref. [1] (called (I) from here 
on), we considered this reaction in the framework of heavy baryon chiral perturbation theory 
(HBCHPT) to third order in the chiral expansion. At that order, the first contributions from 
pion loop graphs, which perturbatively restore unitarity, appear. Besides pion loop diagrams, 
there are tree graphs. Some of these have fixed coefficients, others are accompanied by coupling 
constants not fixed by chiral symmetry. These so-called low energy constants (LECs) must be 
determined by a comparison to data or can be estimated using models. As has been shown in 
refs.[2, 3], these LECs encode the information from higher mass states not present explicitly in 
the effective field theory. There are three important reasons to extend the calculations of (I) 
to fourth order: First, only at that order one has a complete one-loop representation. Second, 
it is known that these fourth order corrections can be large (for a review see ref. [4] and an 
update is given in ref. [5]). Third, only after having obtained an accurate representation of the 
isospin-symmetric amplitude, as done here, one can attack the more subtle problem of isospin 
symmetry violation [6, 7] in low energy pion-nucleon scattering. First steps in the framework 
of HBCHPT have been reported in refs.[8, 9, 10]. 

There have been some interesting new developments since (I) appeared: First, a manifestly 
Lorentz invariant form of baryon chiral perturbation theory was proposed in ref. [11] and some 
implications for the nucleon mass and the scalar form factor to fourth order were worked 
out. The same group is also investigating pion-nucleon scattering in that framework [12]. 
By construction, their approach leads to the correct analytical structure of the pion-nucleon 
scattering amplitude, whereas in the heavy baryon approach special care has to be taken in 
certain regions of the complex plane. Second, a different determination of the dimension two 
and three LECs by fitting the HBCHPT amplitude to the dispersive representation (based on 
the Karlsruhe partial wave analysis) inside the Mandelstam triangle was performed in ref. [13]. 
While that method has the a priori advantage that the chiral expansion is expected to converge 
best in this special region of the Mandelstam plane, it is difficult to work out the theoretical 
uncertainties. Another drawback of this procedure is that only three dimension two LECs 
could be determined with sufficient precision (if one uses the third order HBCHPT amplitude 
as input). This is related to the fact that close to the point v — t — the contribution from 
the third order terms is accompanied by very small kinematical prefactors. 

Here, we will follow the approach used in (I), namely to fit to the phase shifts provided by 
three different partial wave analyses for pion laboratory momenta between 40 and 100 MeV. 
This not only allows for a discussion of the uncertainties due to the input but also gives us the 
possibility to predict the phase shifts at higher and lower momenta, in particular the scattering 
lengths and the range parameters. Furthermore, we are then able to directly compare to the 
third order calculation and draw conclusions on the convergence of the chiral expansion. Of 
course, at the appropriate places we will discuss the relation to the work reported in refs.fll, 13]. 
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The pertinent results of the present investigation can be summarized as follows: 

(i) We have constructed the complete one-loop amplitude for elastic pion-nucleon scattering 
in heavy baryon chiral perturbation theory, including all terms of order q A . It contains 
13 low-energy constants plus one related to fixing the pion-nucleon coupling constant 
through the Goldberger-Treiman discrepancy. Their values can be determined by fitting 
to the two S- and four P-wave partial wave amplitudes for three different sets of available 
pion-nucleon phase shifts in the physical region at low energies (typically in the range of 
40 to 100 MeV pion momentum in the laboratory frame). 

(ii) We have performed two types of fits. In the first one, we fit four combinations of the 
dimension two and four LECs, together with five LECs from the third and five from 
the fourth order. This means that the dimension two LECs are subject to quark mass 
renormalizations from certain fourth order terms. Most fitted LECs are of "natural" size. 
In the second approach, we fix the dimension two LECs as determined from the third 
order fit and determine the corresponding dimension four LEC combinations separately. 
We have studied the convergence of the chiral expansion by comparing the best fits based 
on the second, third and fourth order representation of the scattering amplitudes. The 
fourth order corrections are in general not large, but they improve the description of most 
partial waves. This indicates convergence of the chiral expansion. 

(iii) We can predict the phases at lower and at higher energies, in particular the threshold 
parameters (scattering lengths and effective ranges). The results are not very different 
from the third order study in (I), but the description of the P-waves is improved, in 
particular the scattering length in the delta channel and the energy dependence of the 
small P-waves. The errors on the S-wave scattering lengths are as in (I) since they are due 
to the differences in the partial wave analyses used as input. Our theoretical predictions 
are consistent with recent determinations from pionic hydrogen and deuterium [21]. 

(iv) The pion-nucleon sigma term (at zero momentum transfer) can not be predicted without 
further input since at fourth order a new combination of LECs appears, that is not 
pinned down by the scattering data. Therefore, we have analyzed the sigma term at 
the Cheng-D ashen point. Using a family of sum rules which relate this quantity to 
threshold parameters and known dispersive integrals, we find results consistent with other 
determinations using the various partial wave analyses. 

(v) We do not find any improvement of the chiral description of the so-called subthreshold 
parameters as reported in (I). In some cases, the fourth order prediction is worse than 
the third order one. Since our amplitude is pinned down in the physical region, we do 
not expect the extrapolation to the subthreshold region to be very precise. To circum- 
vent this problem, it is mandatory to combine the chiral representation obtained here 
with dispersion relations, see e.g. ref.[12], or by fitting directly inside the Mandelstam 
triangle [13]. 

The manuscript is organized as follows. In section 2, we briefly discuss the effective La- 
grangian underlying our calculation. All details were given in (I), so here we only spell out 
the new terms at fourth order. Section 3 contains the HBCHPT results for the pion-nucleon 
scattering amplitudes g ± ,h ± to fourth order. The fitting procedure together with the results 
for the phase shifts and threshold parameters are presented in section 4. We also spell out the 
remaining problems related to the sigma term and the subthreshold parameters. The appendix 
contains the analytical expressions for the various threshold parameters. 



2 



2 Effective Lagrangian 

The starting point of our approach is the most general chiral invariant Lagrangian built from pi- 
ons, nucleons and external scalar sources (to account for the explicit chiral symmetry breaking). 
The Goldstone bosons are collected in a 2x2 matrix- valued field U(x) = u 2 (x). We use the 
so-called sigma model parametrisation (gauge). We work in the framework of heavy baryon 
chiral perturbation theory, thus the nucleons are described by structureless non-relativistic 
spin-i particles, denoted by N(x). The effective theory admits a low energy expansion, i.e. the 
corresponding effective Lagrangian can be written as (for more details and references, see e.g. 
[4]) ' ^ 2 2 3 

£cS = C^l + + £wN + £nN + AriV + AriV • • • > (2-1) 

where the ellipsis denotes terms of higher order. For the explicit form of the meson Lagrangian 
and the dimension one, two and three pion-nucleon terms, we refer to (I). The complete fourth 
order Lagrangian is given in ref. [14] and the renormalization is discussed in ref. [15]. For com- 
pleteness, we display here the finite terms from C^ N which contribute to elastic 7rN scattering 

+ (eie + ^39l)v X vW (h x ,h up ) + e 17 [S", S"] [h x „ h\ 

+ (ei 8 - ^ (2 + g\) - ^c 4 ) [S^ S>V [h x „ h vp \ + e 19 ( X+ > (u ■ u) 
11 1 

+e 22 [Dp, [D», ( X+ )]] + e 35 iv»v» (x- V) + ^ («„ X-]) 

+e 37 i[5 , ' t , 5"] [u M , [D v , X-}} + e 38 (x+> (x+> + ens^) + en 6 ( det X + det x f ) 

-i c2 ( ^ ^ > - (12^^ + L ( Ji4 - Jl5 ) )(<"«["• A ^ ) 

"(l2^^ " 8^ ( Jl4 ~ Jl5 ) K^" [W ' A V]) + 32^^' (U ' M [U ' A M) 

+(^ + ^)([^ ttA ]^ +h - c 

- ( 12^ + 32^ C4 - (* + ) ^ + h ' C -) 
+ (l2^ («* " + 2^ 3 ) M** W • ^ + hX -) 

+ (ii + i rf3 )(^^^ + H 

+ (oi^ ~ i C3 ) ^<V*0^ + h.c. 
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" (eib^ + i C2 - i ( Jl4 - Jl5 ) ) • m)Da + h ' c -) 

+ kL^ 1 - i - ) ( [5A ' {u » Kp) Dx + h - c -) 

+ ^ A (^ A l v • A ^1 Da + h - c -) + iL? 9 ^ I s ■ D > K] V ' D+ h - c -) 

~vtrt 9A ( iSV Ku] V ' D + h,c -) " ztrt 9A ( iS * vV [v ' A Kv] V ' D + h ' c ') 

L^SV^ [u • D, V] ^A + h.c 



32m 3 



-(^9A Cl + ^d 16 ) (iS A <x> ' uD x + h.c.) - -^([S* 51 [£>„, (x + >] -Da + h.c.) 
+^ 5 (*LX-, « A ]D A + h.c.) + ^d 18 (S x [v ■ D,x-]D X + h.c.) 

+ ( 16^ + i C4 ) 51 K ' ^ ° v + hX ') " i Cl ^ D " 

—^rgA (iu ■ DS ■ Dv ■ D + h.c.) + — (iS ■ uD 2 v ■ D + h.c.) (2.2) 
m 6 v 7 8m 3 v ' 

5u(^> ■ uS ■ Dv ■ Dv ■ D + h.c.) — - — -gA(iS ■ uv ■ Dv ■ Dv ■ D + h.c.) jiV , 



4m 3 
3 

8m 3 



with v^ the nucleoli's four-velocity, the covariant spin-vector, D^ the chiral covariant 
derivative, = {(d^uu^ + u^d^u), h^ v = [D^u^] + [D v ,u^\ and x± — u^X u ^ ± U X* U en- 
coding the explicit chiral symmetry breaking. Traces in flavor space are denoted by (. . .) and 
X- = X- ~ (X-)/2- We remark that the various parameters like gA,m, ■ ■ ■ appearing in the 
effective Lagrangian have to be taken in the chiral SU(2) limit (m u = rrid = , m s fixed) and 

o o 

should be denoted as gA, m , ■ ■ ■■ Throughout this manuscript, we will not specify this but it 
should be kept in mind. It is also worth mentioning the particular role of the terms ~ en 5) n 6 . 
These operators have no pion matrix elements but are simply contact interactions of the exter- 
nal scalar source with the matter fields and thus contribute to the nucleon mass and the scalar 
form factor. This will be of importance later on. Having constructed the effective pion-nucleon 
Lagrangian to order g 4 , we now turn to use it in order to describe elastic pion-nucleon scat- 
tering. To account for isospin breaking, one has to extend this Lagrangian to include virtual 
photons. This has already been done in [8, 10] and we refer the reader to these papers. For a 
systematic study of isospin violation in the elastic and charge exchange channels, one first has 
to find out to what accuracy the low energy nN phase shifts can be described in the isospin 
symmetric framework.* 5 This is the question which will be addressed in the remaining sections 
of this investigation. 

#5 An investigation of isospin violation in the threshold amplitudes to third order was reported in ref.[9]. 
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3 Pion— nucleon scattering 



3.1 Basic definitions 



In this section, we only give a few basic definitions pertinent to elastic pion-nucleon scattering. 
For a more detailed discussion, we refer to (I). In the center-of-mass system (cms), the amplitude 
for the process n a (qi) +N(p 1 ) — > rr b (q 2 ) + N(p 2 ) takes the following form (in the isospin basis): 



rpba 
1 ttN 



E + m 
2m 



+i e bac r c 



g + (uj, t)+ia- (q 2 x q\ ) h + (u, t) 
g~(uj,t) + ia ■ (q 2 x q x ) h~(u,t) 



(3.1) 



with u — v ■ qi 

• 2 - a 2 ~ 

q 2 



E = (<f 



m 



I 1 / 2 the nucleon energy 



v ■ q 2 the pion cms energy, E\ = E 2 
and q \ = q 2 = q L = ((s — M 2 — m 2 ) 2 — Am 2 M 2 ) / '(4s) . t = {qi—q2) 2 is the invariant momentum 
transfer squared and s denotes the total cms energy squared. Furthermore, g ± (u,t) refers to 
the isoscalar/isovector non-spin-flip amplitude and /i ± (u;,t) to the isoscalar/isovector spin-flip 
amplitude. This form is most suitable for the HBCHPT calculation since it is already defined 
in a two-component framework. 



3.2 Chiral expansion of the amplitudes 

What we are after is the chiral expansion of the various amplitudes g ± ,h ± . These consist of 
essentially three pieces, which are the tree and counterterm parts of polynomial type as well as 
the unitarity corrections due to the pion loops. To be precise, we have 

X = X trcc + X ct + X loop , X = g ± ,h ± , (3.2) 

where the tree contribution subsumes all Born terms with fixed coefficients, the counterterm 
amplitude the ones proportional to the dimension two, three and four LECs. The last term in 
Eq.(3.2) is the complete one-loop amplitude consisting of terms of order q 3 and q 4 . The latter 
is a complex-valued function and restores unitarity in the perturbative sense. Its various terms 
are all proportional to 1/E 4 . We remark that the topologies of the new loop graphs are not 
different from the ones already present at third order. The loops of fourth order have exactly 
one insertion from the dimension two Lagrangian. Note that we treat the chiral symmetry 
breaking scale A x ~ 1 GeV on the same footing as the nucleon mass. In principle, one could 
also organize the loop expansion, which proceeds in powers of 1/-F 2 , and the 1/m expansion 
independently of each other (with some prescription for the mixed terms). These amplitudes 
are functions of two kinematical variables, which we choose to be the pion energy and the 
invariant momentum transfer squared, i.e. X = X(u,t). In what follows, we mostly suppress 
these arguments. The full one-loop amplitude to order g 4 is obtained after mass and coupling 
constant renormalization, 

(g A , rh, F, M) -> (g A , m, F w , M w ) . (3.3) 
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To fourth order, these read (we also give the corresponding Z-factors for the pion and the 
nucleon) 

2M 2 A 

z i 9 2 a sm 2 g 2 A 9M 3 f35 ^ 

Ml = M 2 {l + ^3 + ^) , (3.6) 



with 



and 



m = m -AM ci 



F 2 •> 2F 2 
3g 2 A M 3 



32vrF 2 



1 3M 4 c 2 Sg 2 A M 4 . . 

— M (16e 38 + 2e 115 + -e 116 ) + - , (3.7) 

= F|l + ^ 4 -^} , (3.8) 

^ ^/i M % + 4M2 r/ m+ ^ (a m ^ 

^ " Fl 1_ ^ 4 + ^7 dl6(A) + 4^1, A7r "4^J 

-^T- + c 3" 2c 4 + — --^ , (3.9) 



67rF 2 y 8m " V ! 2m) Am 

1 . Mir 



A n = 2Ml(L + — \n^)+0(d-A) (3.10) 

L = ^(^ + I(W-ln4.)), (3.11) 

where Euler's constant 7# = 0.557215 has been used, A is the scale of dimensional regularization 
and d the number of space-time dimensions. 

After these preliminaries, we give the final expressions for the tree, counterterm and loop graphs 
at fourth order in terms of the renormalized quantities: 

Tree and counterterm graphs: 



m 3 F^g+(u,t) = 



9 2 A 



64w 4 



32c/M 2 t + 45M 4 t 2 - llM 2 t 3 - 80M^ + 52M* 

+7u 2 t 3 + llOcu 2 M 4 t - A9uj 2 Mlt 2 - 7Qu 2 Ml + lluA 
+t 4 - AujH + 28c/M 41 



+\mI{2u 2 - 2Ml + t)m Cl + 8M 4 m 3 Cl || 

2 7T 

i(-22c; 2 M 2 + 8M 4 + SuH + 14c; 4 - 4M 2 t)mc 2 



+ 



+ i(-4cu 2 M 2 + 2uH + AMI ~ 4Mlt + t 2 )mc 3 - 16M 2 cj 2 m 2 Cl c 2 
8 



m 3 F 2 h + (uj,t) 



m 3 F 2 g-(u,t) = 



m 3 F*h-(LU,t) 



+uj 2 tm 2 (d u - d 15 ) + ^M 2 (4M 4 + t 2 + AuH - AMlt)m 2 d l8 

-4(4M 2 t - 4M 4 + t 2 )m 3 e 14 - 8u 2 (-2M 2 + t)m 3 e 15 
+ 16c; 4 m 3 e 16 - l.U:( 2M: + t)m 3 (2e l9 - e 22 - e 36 ) 



+16a> 2 M 2 m 3 (e 20 + e 35 - 



9Ad 



16 - 



8m 



) + 8M> J (e 22 - 4e 38 ) , 



9\ 



32c; 4 



5c/r + 27M*t - 9M z t z - 25u z M;t - 28M^ + t d 



' u - 9M 2 t 2 
-m.^ Ml - 4c; 4 M 2 + 3ujH 



; 2 M , 
1 



1 



+M;mci - -u 2 mc 2 + -(-2M; + t)mc 3 
+i(8c; 2 - 4M 2 + t)m 2 (J 14 - J 1B ) - 0M 2 (-4M 2 + t)m 2 rf 18 



64c; 4 



32o> 4 M 2 t + 45M 4 t 2 



llM 2 t 3 - 82M^t + 56M* 



+7u/r + 110c; 2 AO - 49c/M 2 t 2 



4.2 



+8o/ + F- 2J°t - 8c; b M 2 + 24c; 4 M 4 



1 



4c; 8 + 8c; 6 M 2 - 4c; 4 M 4 - uoH + c; 4 M 2 t 
+^t(4cu 2 - 4M 2 + t)mc A 

-i(-8c; 2 M 2 + 4c; 2 t + 8M 4 + t 2 - 6M 2 t)m 2 (cJ 1 + J 2 ) 
+3c; 2 (4c; 2 - 4M 2 + t)m 2 rf 3 + 2M 2 (4c; 2 - 4M 2 + t)m 2 d 5 



9a 
"4c; 2 
~9 2 A 



32c; 4 



M 2 (t 2 + 2uo 2 t - 8a; 4 + 8M 4 - QM 2 t)m 2 d 18 , 

5tuH 2 + 27M 4 t - 9M 2 t 2 - 25u 2 M 2 t - 26M 6 + t 3 

u\uo 2 - Ml) 



+30a; 2 M 4 - 4c/M 2 + 3c/t - 2c; 6 



16c; 4 



(3.12) 



(3.13) 



(3.14) 



-M 2 m Cl + -(2c; 2 - 2M 2 + t)mc 4 + ^-M 2 {2u 2 - 2M 2 + t)m 2 d 18 
8 2ui 



4(-2M 2 + t)m 3 e 17 + 8u 2 m d e 18 + 4M 2 m 3 (2e 21 - e 37 ) . 



(3.15) 



Loop graphs: 
mF*g + (u,t) 



12c; 3 



JoH 2M 2 g\(tM 2 - 2c; 4 + 4c; 2 M 2 - too 1 - 2M 4 ) 



12c; 3 



+cu 2 g 2 A (-12c; 2 M 2 - M 2 t + c; 2 t + 12c; 4 ) + 6c; 4 M 2 

^J (-cj) [2g\{<6Ml + 8c; 6 - 10c; 4 M 2 + 3tc; 2 M 2 + M 2 t 2 + 2tc; 4 

-4c; 2 M 4 - 5tM 4 ) - uj 2 g 2 A {M 2 t + 4c; 2 M 2 + c; 2 t - 4c; 4 ) 
+6c; 4 (-3M 2 + t + 4c; 2 )" 



7 



mF*h + {uj,t) = 



mF^g (u,t) = 



1 dJ 
12c; 2 Ouj 



u) g\{8^ 



u 2 t 2 - 8M° - M 2 t + 2Au 2 M* + 6tc; 4 + 6tM* 



-12tco 2 M 2 



2Au A Ml) 



+2u 2 g 2 A (Acu 4 - Mlt - Suj 2 MI + cuh + AM*) 
+3c/(4c; 2 -AMl + t) 

+§!^(0,t)(12A# - AMI ~ $M 2 t 2 + 2t 3 ) 
32 oui 

~Ht)[39At(2t - Ml) - 48M> Cl (M 2 - 2t) 

+2mc 2 (2t 2 + AMI ~ 9*M*) + 12mc 3 (-5tM 2 + 2t 2 + 2M 4 

1 



{ , j\2g\(-r>2j'M; + 24c; 7 + uH 2 + lOtc; 5 - 96ttM 7 + 96ttM> 2 



r2,2 



+72nM°t - 12ttM^ - 48ttM>^ - 12wM*t + 6wM# 
+ llc; 3 M 2 t + 46M 4 c; 3 ) 



+g 2 A (-18coH 2 - 288c; 5 M 2 + 60c; 3 M 4 + 33c; 3 tM 2 + 192c; 7 
+64cA) + 4c; 3 mc 2 (2t 2 + 6M 4 - 13M 2 t) 



(3.16) 



12c; 3 



J (^)(cu 2 - M 2 ) [<? 2 M 2 - 2c; 2 + 8mu 2 {c 3 - c 4 ) 



^sM-w) [9 2 a(M 2 uj 2 + Mlt - 3M 4 + 2c; 4 ) 

+ (lu 2 - M 2 )(2c; 2 - 8c; 2 m(c 3 - c 4 )) 

+jk^(-oo)(Au A - M 2 n t - 8u 2 Ml + cA + 4M 4 ) 



24c; 2 dcu 

g\dK 



32 a7 (0 ' *)(" 9J ^* + 4M * + 2t2) " Y m{2t ~ M " ] 



9a 



1152w% 2 



2g A (-12irM* - 6ttM> 2 + 6vrM;t + 32c; + 4c/t 

-16c; 3 M 2 + 12c;M 4 - 3c;tM 2 ) 
+3c; 2 (3c;M 2 - 6ut - 16vrM 3 ) + 192tt M 3 c; 2 m(c 3 - c 4 ) 



(3.17) 



24c; 3 



24c; 3 



JoH [<?1(^ 2 M 2 + 2c; 4 M 2 - 4c; 2 M 4 - tM 4 + 2M 6 ) 

+g 2 A (-Quj 2 M* + 12c; 4 M 2 + 3c; 2 M 2 t - 3c/t - 12c; 6 

+8c; 2 m(c 3 - c 4 )(2c; 4 - 4M 2 c; 2 + 2M 4 - tM 2 + c; 2 t)) 
+6c; 4 (-M 2 - 16M 2 mci + 8c; 2 m(c 2 + c 3 )) 

Jo(-oj) U(-5tM 4 - 10c; 4 M 2 - 4c; 2 M 4 + 3tc; 2 M 2 + 6M 6 



+8c; 6 + Mlt 2 + 2tc; 4 ) 
+^c; 2 (6M 4 - 4c; 2 M 2 + c; 2 t - 3M 2 t + 4c; 4 

+8m(c 3 - c 4 )(-2M 4 + tMl - uo 2 t - 2c; 4 + 4c; 2 M 2 )) 
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1 dJ 

48c; 2 dco 



+6c; 4 (-3M 2 + t + 4c; 2 + IQMlmd - 8u 2 m(c 2 + c 3 )) 



{-u) [g\(-\2tuj 2 Ml - 24c; 4 M 2 + 6tM 4 + 8c; 6 + 24c; 2 M 4 

-8M 6 + 6tc; 4 - M 2 t 2 + oo 2 t 2 ) 
+4^ (4c; 4 + 4M 4 + ojH - 8u 2 M 2 - M 2 t) 

+6c; 4 (-4M 2 + 4c; 2 + *)" 



9\ 



+^ (0,tV(4M 4 + t 2 -4M 2 t) 

o 



4^(0,t)(-10M 2 t 2 + 32c; 2 M 4 - 24c; 2 M 2 t + 4c; 2 t 2 + 32M 4 t 



64 du 



-32M 6 + 1 3 ) 



+^h(t) [2g\{t 2 - 4M 2 c; 2 + 4M 4 + 4cA - 5M 2 t) 
-8M 2 t + 16M 4 + t 2 - 16c; 2 M 2 + 4c; 2 t 

1 



2304c; 3 7r 2 



48ttM^ - 12ttM^ 2 
36c; M 4 t + 6c;M 2 t 2 



c/l(-487rMj - 12ttM 3 cA - 24ttM 5 c; 2 

+96c/ + AujH 2 + 40tc; 5 + 72ttc; 4 M 3 
-12c; 3 M 2 t + 32M> 3 + 48c;M 6 - 176c; 5 M 2 ) 
+#^(7cA 2 + 72vrM 3 cA - 144ttM 5 c; 2 + 72nuj 4 tM 7T 
-456c/M 2 + 264c; 3 M 4 - 94c; 3 tM 2 + 192c/ + 76c; 5 t 



+216ttM 3 c; 4 + 1927rM 3 c; 2 m(c3 - c A ){-t - 2c; 2 + 2M 2 )) 



r3, ,2 



9\ 



12c; 3 



+2c; 3 (-10tM 2 - 24c; 2 M 2 + t 2 + 24M 4 + 4c; 2 t) 
J Q {uo){uo 2 - M 2 )(-2M 2 ^i + 3c; 2 + 8c; 2 mc 4 ) 



(3.18) 



9\ 



+J^3 J o(-u) [2^i(-3M 4 + 2c; 4 + M 2 t + M 2 oo 2 ) 



9a dJ 
12c; 2 dcu 



-uo 



9\ 



~I (t) [2si(-5M 2 + 2t) + (AM 2 -t)(l + Amc A ) 



WV 2 - 3M 2 + 8mc 4 (u 2 - Ml)) 

g 2 A {-8Mlu 2 - Mlt + ooH + 4c; 4 + 4M 4 ) 

+2c; 2 (c; 2 — M 2 )" 

9 2 A dK 



K (0,t)uj(-4M 2 + t) + 



32 du 



(0,t)(8M 4 -6M 2 t + t 2 ) 



24 



1 



1152c; 3 7r 2 



[g\(96nM^ - 64c; 3 M 2 + 12c;tM 2 + 8cA - 24vrM 3 t 

-48vrM 3 c; 2 + 64c; 5 ) 
+^(-18c; 3 M 2 + 48c; 5 + 36nu*M w + llc; 3 t 
+ 16c; 3 mc 4 (15M 2 - 4c; 2 )) 

+2c; 3 (6M 2 - * + 4mc 4 (6M 2 - f)) 



(3.19) 
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We have used the loop functions of ref.[16]. The are scale-independent (using the same 
procedure to eliminate the chiral logarithms as detailed in (I) for the di). It is important to 
stress that we have obtained a more precise representation of the imaginary parts as compared 
to (I) since in that paper, only the leading terms were included. Here, we also have the next- 
to-leading order corrections of the unitarity corrections. Clearly, unitarity is perturbatively 
fulfilled, i.e. Im ~ (Re T^) 2 in a highly symbolic notation, where refers to the chiral 
representation of the 7rN amplitude to n th order. It goes without saying that we also expect 
the corresponding real parts to be given more accurately. 



3.3 Counterterm combinations 

In the counterterm and loop contributions given above, a set of LECs appears. Most of these 
only enter in certain combinations and some only lead to quark mass renormalizations of the 
dimension two LECs. Therefore, it is instructive to work out how many independent local 
contact terms can contribute to 7rN scattering to fourth order. This can be most easily done 
based on a dispersive analysis by counting the number of possible subtractions. For that, 
it is most appropriate to describe the pertinent T-matrix in terms of the standard invariant 
amplitudes A and B, 

T± N = A ± +</B ± , (3.20) 

in a highly symbolic notation. The invariant amplitudes are functions of two variables, which 
one can choose to be v and t; these count as 0(q) and 0(q 2 ), respectively. The most general 
polynom for the four amplitudes A ± ,B ± to fourth order commensurate with crossing and the 
other symmetries thus takes the form 

Ap ol = a~l + a\t + a^v 2 + a\ t 2 + a^tis 2 + z/ 4 , 

B pol = b- 1+ b 2 t + b^v 2 . (3.21) 

Certain combinations of dimension two, three and four LECs are related to the subtraction 
constants (af, . . . , 63 ). We refrain from giving the precise relationship here since it is not 
needed in what follows. Therefore, in total we have 14 LECs since at third order there is one 
more related to the Goldberger-Treiman discrepancy, i.e. a local term with a LEC which allows 
to express the axial-vector coupling g A in terms of the pion-nucleon coupling g nN , i.e. 

_g A m( 2M 2 d 18 \ 
9^-^-{l —). (3.22) 

This term is important if one wants to properly account for the Born terms expressed as a 
function of the pion-nucleon coupling constant. If one then calculates to orders q 2 , q 3 and q 4 , 
one has to pin down 4, 9 and 14 LECs, respectively. This pattern is quite different from the total 
number of terms in the Lagrangian allowed at the various orders (7, 23, and 118, respectively); 
it is a general rule that simple processes do not involve an exorbitant number of LECs. Indeed, 
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the terms proportional to the dimension four LECs (i = 19,20,21,22,35,36,37,38) only 
amount to quark mass renormalizations of the dimension two LECs q (i — 1,2, 3, 4) via 



ci = ci -2M 2 (e 22 -4e 38 ) , 



c 2 = c 2 + 8M 2 (e 20 + e 35 - 




) 



c 3 = c 3 + 4M 2 (2ei 9 - e 22 - e 36 ) , 
c 4 = c 4 + 4M 2 (2e 2 i - e 37 ) . 



(3.23) 



We have also used these parameters in the one-loop graphs of order q 4 , although this induces 
some higher order contributions. This is a very general phenomenon of CHPT calculations at 
higher orders (for a discussion, see e.g. ref.[17]). There are different ways of determining the 
LECs. As in (I), we use data from the physical region for doing so. Our first strategy is to fit 
the renormalized q, called q here, together with the four (combinations of) dimension three 
LECs di + J 2 , ds, d§, du — di$ and di$ and the genuine dimension four LECs &u, £15, ei6, ei7, eis- 
As enumerated before, we thus have 14 free parameters. In such a fit, we cannot disentangle 
the Cj into their quark mass dependent and independent pieces without further information 
from other processes. This defines our best fit. To study the convergence compared to the 
lower order calculations, we also show the best fit from (I) and a best second order fit based 
on tree diagrams with the dimension two insertions ~ q. One can argue that the second order 
contribution is given by the amplitude up to second order, with the c's taking on their values 
as given by the best fit at that order. By including the amplitude at third order, the values of 
the c's will change, but these changes are considered to be effects of third order. (The same is 
valid of course for the c'd and gPs, when going from third to fourth order.) One might also be 
interested in how big the contribution from genuine second and third order terms are (terms 
really proportioanl to q 2 , respectively to q 3 ). In order to address this question, we consider an 
alternative strategy, in which we fix the q (i = 1,2,3,4) to their values determined from the 
best fit up to third order and use the four combinations of dimension four LECs appearing in 
eq.(3.23) as fit parameters. Of course, this leads to the same number of free parameters, but 
this second method allows for a clean separation of the contributions from the various orders. 
Clearly, physical observables do not depend on this reshuffling of fit parameters (modulo higher 
order corrections effectively included when using the q in the loop graphs). 

4 Results 

4.1 The fitting procedure 

There are various possibilities to fix the LECs, a general discussion is given in (I). We proceed 
here along the same lines as in (I), namely we fit to the phase shifts given by three different 
partial wave analyses in the low energy region. As input we use the phase shifts of the Karlsruhe 
(KA85) group [18], from the analysis of Matsinos [19] (EM98) and the solution called SP98 
from the VPI/GW group [20]. * 6 In contrast to what was done in (I), we do not assign a 
common error of 3% to the Karlsruhe and VPI phases, but rather mimic the uncertainties 

# 6 In the meantime, novel solutions like SM99 have appeared. Since these are not very different from SP98 
and we want to have a direct comparison with the results of (I), we use SP98 here. We come back to this later. 
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of the Matsinos analysis in all cases, which is 1.5% for S31, 0.5% for Su, 1% for P 33 and 
3.5% for the other P-waves. This assignment gives more weight to the better determined 
larger partial waves and is more natural than one common global error.* 7 The LEC dig 
is fixed by means of the Goldberger-Treiman discrepancy, i.e. by the value for the pion- 
nucleon coupling constant extracted in the various analyses. The actual values of g n N are 
g nN = 13.4±0.1 , 13.18±0.12 , 13.13±0.03 , for KA85, EM98 and SP98. Throughout, we use 
g A = 1.26, F w = 92.4 MeV, m = 938.27 MeV and M n = 139.57 MeV. Finally, we remark that we 
do not use the value of the pion-nucleon a-term in the fitting procedure. This has two reasons: 
First, as noted before, we only want to use information from the physical region to pin down 
the LECs and second, it is known that the convergence of the chiral series for this quantity is 
slow [17]. Before presenting the results of the actual fits, we already anticipate that the EM98 
data basis will lead to the smallest x 2 f° r the following reasons. First, this data base is only 
covering the low-energy region of pion-nucleon scattering. Also, the representation is available 
on a denser grid of points in momentum transfer. In contrast, the KA85 and SP98 analyses 
span a much larger range of energy and thus uncertainties also from higher energies will play 
a role in the energy range considered here. Furthermore, in (I) we had already discussed that 
the extraction of the threshold parameters from the SP98 analysis is not unproblematic. Note, 
however, that the model underlying the EM98 analysis should not be used in the unphysical 
region, quite in contrast to the dispersion theoretical approach on which the KA85 phase shifts 
are based. 

4.2 Phase shifts and threshold parameters 

After the remarks of the preceding paragraph, we can now present results. For the KA85 case, 
we have fitted the data up to 100 MeV pion lab momentum (i.e. 4 points per partial wave at 
q n = 40, 60, 79, 97 MeV). For the analysis of Matsinos, we use 17 points for each partial wave in 
the range of q n = 41.4 — 96.3 MeV. For the VPI SP98 solution, we use the 5 data points in the 
range between 60 and 100 MeV, which give a stable fit. Of course, we could now extend the 
fits to higher energies than it was done in (I), but for a better comparison we do not show the 
results of these extended fits here. As discussed in section 3.3, we have two options for pinning 
down the LECs. Using strategy one, i.e. working with the q, we call the fits corresponding 
to the Karlsruhe, Matsinos and VPI analysis, fit 1, 2 and 3, in order. The resulting LECs are 
given in table 1. Note that the error on the LECs is purely the one given by the fitting routine 
and is certainly underestimated. We remark that the q and di (or combinations thereof) are 
mostly of natural size, whereas some of the come out fairly large. Also, there is some sizeable 
variation in the actual values of most LECs among the different fits. The resulting S- and P- 
wave phase shifts are shown in figs. 1 (fit 1), 2 (fit 2) and 3 (fit 3), in order. The corresponding 
X 2 /dof is 0.40, 0.008 and 0.44 for fits 1,2 and 3, respectively. The description of the phase shifts 
is excellent for fit 2. For fits 1 and 3, the S31 and P\\ partial waves at pion momenta above 
150 MeV are somewhat off. In all cases, the description of the P33 partial wave is improved as 
compared to the third order calculation. Since we do not fit data below q n = 40, 41, 60 MeV 
(fit 1,2,3) , the threshold parameters are now predictions. These are shown for the various fits 

course, this might induce some mismatch in the sense that real errors associated to the KA and VPI/GW 
phases are different from the ones of the Matsinos analysis. We believe, however, that this procedure is preferable 
to the one using common global errors. 
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Table 1: Values of the LECs in GeV 1 , GeV 2 and GeV 3 for the q, di and e^, respectively, 
for the various fits based on the first procedure as described in the text. 

in table 2, in comparison to the empirical values of the various analyses. First, we observe that 
the numbers resulting from the one-loop calculation are consistent with the "empirical" ones. 



Obs. 


Fit 1 


Fit 2 


Fit 3 
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0.45 
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7.71 
1.78 
-5.87 
13.04 
-1.17 
-8.21 


0.27 
-7.16 
8.67 
1.15 
-4.81 
13.18 
-0.79 
-7.54 
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-4.40 
9.17 
0.77 
-5.53 
13.27 
-1.13 
-8.13 


0.41 ±0.09 

-4.46 
7.73 ±0.06 

1.56 
-5.46 ±0.10 
13.13 ±0.13 
-1.19 ±0.08 
-8.22 ±0.07 


0.0 ±0.1 
-4.83 ±0.10 
8.83 ±0.07 
0.07 ±0.07 
-5.33 ±0.17 

13.6 ±0.1 
-1.00 ±0.10 
-7.47 ±0.13 



Table 2: Values of the S- and P-wave threshold parameters for the various fits as described 
in the text in comparison to the respective data. Note that we have extracted 6q+ from the 
Matsinos phase shifts and thus no uncertainty is given. Units are appropriate inverse powers 
of the pion mass times 10~ 2 . 

The agreement with the threshold parameters based on the chiral amplitude with the ones 
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based on the approaches underlying the various partial wave analyses is best for fit 2, slightly 
worse for fit 1 and clearly problematic for some of the parameters of fit 3. The reason for this 
was already spelled out in (I). The bands for the S-wave scattering lengths a^~ + and ag + are as 
in (I) since the uncertainties extracted there are mostly due to the input and not to the theory. 
They agree with the recent determinations from the shift and width of pionic hydrogen and 
deuterium, cf. Fig. 2 in ref.[21]. For comparison, we translate our bands on the isoscalar and 
isovector scattering lengths into the physical ones, 

a^-p^on = -0.131... -0.117 M~ l [(-0.128 ± 0.006) M~ l ] , 

a^-p^-p = 0.073... 0.093 M' 1 [(0.0883 ± 0.0008) M~ l ] , (4.1) 

where the experimental numbers (in the square brackets) are taken from ref.[21]. Note, however, 
that recent progress in calculating n~p atoms in effective field theory lets one expect that the 
uncertainty due to electromagnetic corrections for the band derived from the hydrogen shift has 
been underestimated, see e.g. ref. [23]. Furthermore, only recently deuteron wave functions have 
been obtained precisely enough in an EFT approach to readdress the question of the deuterium 
shift constraining the elementary nN amplitudes. It would also be worthwhile to repeat the 
EFT calculation of pion-deuteron scattering [22] using our fourth order nN amplitudes as 
input. As argued before, we can study the convergence of the chiral expansion. In figs. 1-3, 
the dot-dashed, dotted, dashed and solid lines refer to the best fits up to first, second, third 
and fourth order, respectively. Since we used the errors of the Matsinos analysis, it is best to 
consider fit 2 shown in fig. 2. In most cases, the fourth order corrections are smaller than the 
third order ones, indicating convergence. This could not be concluded from the third order 
calculation, compare the discussion in (I) and ref. [16]. Note also that in some partial waves 
the second order result is close to the data. The resulting values of the q are very different 
from the ones given in table 4. The second order best fit based on the (KA85, EM98, SP98) 
analysis leads to c x = (-0.81, -0.77, -1.06), c 2 = (2.47, 2.69, 2.36), c 3 = (-3.78, -3.96, -4.04) 
and c 4 = (2.49,2.64,2.35) (all in GeV" 1 ).* 8 That these values are very different from the 
ones based on a one-loop third order amplitude fit was already pointed out in ref. [24]. It is of 
particular interest to study the convergence of the S-wave scattering lengths, which has been 
already discussed in ref. [25, 26] estimating LECs from resonance saturation. Our results are 
summarized in table 3. Although it was already shown in ref. [26] that there are no fourth order 
corrections to aQ + , the readjustment of the LECs when going from third to fourth order leads 
to a small difference. That this difference is so small is also in agreement with ref. [26] , where it 
was argued that the dominant correction to the Weinberg-Tomozawa low-energy theorem is a 
pion loop effect. For fits 1 and 2, the fourth order correction to the isoscalar S-wave scattering 
length is fairly small, even for fit 3 the dominant correction is the one from second to third 
order. 

The second option is to keep the dimension two LECs fixed to their value determined from 
the third order fits and fit the additional four dimension four combinations. This allows for a 
clean discussion of the various contributions to the chiral expansion. The results for the LECs 
are shown in table 4.* 9 The stability of the values for most of the LECs is better than in the 
previous case and the fourth order LECs e^, i — 14, . . . , 18 smaller (more natural). Some of the 

* 8 Note that for such a second order fit the differentiation between the Ci and becomes meaningless. 
#9 The slight differences for the values of the as compared to the ones given in (I) stem from the fact that 
we use different error bars for the KA85 and SP98 partial waves as explained before. 
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Table 3: Convergence of the S-wave scattering lengths for the three fits. 0(q n ) means that all 
terms up-to-and-including order n are given. Units are 10~ 2 /M n . 



additional combinations of dimension four LECs are fairly large and vary considerably for the 
various fits. The resulting S- and P-wave phase shifts are shown in figs. 4 (fit 1*), 5 (fit 2*) 
and 6 (fit 3*), in order. The corresponding x 2 /dof is 0.50, 0.14 and 0.58 for fits 1*,2* and 
3*, respectively. In these plots a different way of looking at the convergence properties of the 
amplitude is adopted: all four curves are based on the same fit and are thus obtained with 
the same set of LECs. The dot-dashed, dotted, dashed and solid lines show the contributions 
from the amplitude up to first, second, third and fourth order, respectively. We remark that 
the fourth order contributions are mostly small, with the exception of the Pu and P 13 partial 
waves. The threshold parameters determined from these fits come out very close to the values 
given before and we thus refrain from adding another table here. Similar remarks hold for the 
convergence of the S-wave scattering lengths, only that in this way of fixing the LECs there is 
indeed no contribution to the isovector one from fourth order. 

4.3 Sigma term, subthreshold parameters and further comments 

In this section, we briefly discuss the so-called subthreshold parameters and the sigma term. 
Already in (I) we noted that the representation of the chiral amplitude, when pinned down by 
scattering data, is not very precise in the unphysical region. In particular, the small isoscalar 
amplitudes are obtained from various contributions, which are individually much larger than 
their sum. Consequently, this fine balance which is enforced through the fit in the physical 
region down to the scattering lengths is disturbed because the strict 1/m expansion performed 
here does not properly account for all cuts appearing in the 7rN amplitude. In fact, using 
our fourth order representation, we do not find an improvement of the subthreshold parame- 
ters as given in (I), in some cases even a clear disimprovement. This problem could e.g. be 
circumvented in the formulation of ref.fll]. Another option is to pin down the LECs inside 
the Mandelstam triangle [13], which will lead to an improved representation in the unphysical 
region. This is also reflected in the prediction for the sigma term, which came out rather large 
in the fits shown in (I), but was considerably different (and consistent with the result from 
dispersion theory) based on the method used in [13]. 

We now consider the sigma term, which is the matrix element of the explicit chiral symmetry 
breaking part of the QCD Lagrangian sandwiched between proton states at zero momentum 
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Table 4: Values of the LECs in GeV -1 , GeV -2 and GeV~ 3 for the q, di and e,, respectively, 
for the various fits based on the second procedure as described in the text. The * denotes an 
input quantity. 



transfer. While at third order we can directly give the sigma term, cr(0), in terms of the LEC 
Ci, this can no longer be done at fourth order due to the appearance of the LEC combination 
2eii 5 ±e 116 /2. These operators contribute to the nucleon mass shift and the sigma term (scalar 
form factor) as noted before. These contact interactions have no pion matrix-elements and 
therefore can not appear in the scattering amplitude, even not in higher order loop graphs. 
We therefore use a more indirect method to determine the sigma term. For that, we consider 
£ = F^D + {v = 0, t = 2M 2 ) which can be related to <j(0) by the venerable low-energy theorem 
of ref.[27]. There exists a whole family of relations between £ and certain combinations of 
threshold parameters, as detailed in ref. [28]. These relations have been worked out to third 
order and should be generalized to fourth order. We will use here the version given in ref. [29], 

£ = ttF 2 [(4 + 2/i + n 2 )at + - 4M 2 5+ + + 12//MX+] + , (4.2) 

with £ — ~ 12. 6 Me V and // = M n /m ~ 1/7. Using the pertinent threshold parameters from 
the fourth (third) order representation, we find £ = 65 (62) MeV, 73 (79) MeV and 90 (82) MeV 
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for fits 1, 2 and 3, respectively. A special variant, which also contains some fourth order pieces, 
has recently been given by Olsson [30], 

E = [F*F(2M 2 J], (4.3) 

F(2Ml) = 14.5(4+ - 5 - 06 ( a l+ ) 2 ~ 10 - 13 («o+) 2 - 16 - 65 K+ ~ °- 06 at- + 5 - 70 a t+ ~ 0-05 , 

with the quantities on the right-hand-side being given in units of the pion mass. This leads to 
E = 73 (69) MeV, 85 (91) MeV and 104 (93) MeV for fits 1, 2 and 3, respectively. We consider 
the differences between the results based on eqs.(4.2) and (4.3) (and also using the fourth 
order results for the threshold parameters in the third order representation, eq.(4.2)) as an 
indication of the size of the fourth order terms. We note that the values we find for the 
Karlsruhe analysis are consistent with the direct determination based on hyperbolic dispersion 
relations [31], whereas the results based on the SP98 partial waves lead to a sizeably larger 
value than advocated by the VPI/GW group [32]. 
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A Threshold parameters 

In this appendix, we give the analytical expressions for the threshold parameters up to fourth 
order. These read: 



Affi-fl% + 8m(-2ci + C2 + c 3 )] 
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+ 167rm(m + M 7r )F2 + 30727r 2 m(m + M„)F* 
g 2 M 2 (-16m 2 - 2mM w + M 2 ) M 2 c 2 

2567r(m + M w )m 4 F 2 + 2vr(m + M,)mF x 2 
M 2 Cl c 2 (4m 2 - 2mM T + M 2 ) _ M^c 1 £ 3 (2m - mJ) 
rr(m + M^)m 2 F 2 " 2n(m + M n )mF* 

M 2 [8m 2 (J 14 - d 15 ) - g A dis(M 2 - 2mM 7T + Am 2 )] 

16vr(m + M w )m 2 F 2 
M 2 (e u + ei 5 + e 16 ) (-8m 2 + 2mM 7T - M 2 ) M 2 (2mM 7r - M 2 ) (e 22 - 4e 38 ) 
n(m + M v )mF 2 27r(m + M 7r )mF 2 

M 2 [2(e 20 + e 35 - ^) + (2e 19 - e 22 - e 36 )](-4m 2 + 2mM 7T - M 2 ) 

2-(;// • .1/, 

M 2 

92167r J (m + M n )m z F£ 
+g\(m 2 (-88 + 192tt) - 36mM^ + 18M 2 ) + mci(1248m 2 - 144mM 7r + 72M 2 ) 
-24m 3 c 2 + mc 3 (-768m 2 + 72mM n - 36M 2 )] , (A.3) 
2m 2 - 2mM n + M 2 
32irmM n (m + M n )F 2 



lQn(m + M W )F 2 

M„g A (-10m? - 2mM n + M 2 ) M^c 4 

1287rm 3 (m + M^)F2 4vr(m + M n )F 2 

f M^r__ 

47rm(m + M W )F 2 



izs7rm^m -+- ivi n )r^ wrym -\- ivi n )r£ 

M v [{di + d 2 + d 3 )(6M 2 - 2mM T + M 2 ) + J 5 (4m 2 - 4mM 7r + 2M 2 )] 

47rm(m + AOF 2 
M 7r (4m 2 + QmM w - 3M 2 - Ug A m 2 ) 
7687r 3 m(m + M n )F£ 
3g\M 2 M 2 (d l + d 2 + 3d 3 + 2d 5 + g A d 18 ) 

647r(m + M^)m 2 F 2 * ' 2n{m + M W )F 2 

M 2 (18 + 69irg 2 A +(4-12n)g A ) 

23 047r 3 (m + M^F* ' 1 ' ' 



g\m 



~24irM 7r (m + M W )FZ 

3-2g A + 8mc 4 
~487r{m + M n )F 2 
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~" 967r(m + M n )F 2 
M^m 2 [3 + g\(21 + 36tt) + g\{2 + 24tt)] 
34567r 3 (m + M V )F* 
9a m I Mid , M*(d 1 + d 2 + 3d 3 + 2d 5 + g A d 18 ) 



1927r(m + M n )m 2 F% Qn(m + M n )mF% Qn(m + M n )F% 

2mMl(2e 17 + 2e ls + 2e 2 i - e 37 ) 



37r(m + M 7r )F2 
M 2 

- [18 + g\(12 - 33tt) + g\(A + 30tt) + 967r^rac 3 



Q9l2n 3 (m + M W )F* [ 

•(176 - 967r)^mc 4 ] , (A.5) 



mM^m + M^F 2 
#i + 2mc 3 



127r(m + A^)F 2 

M^[c 2 + 2m(d 14 - d 15 + g A di 8 )} mM 7r ^(23l7r + ^(112 + 96tt)) 



67r(m + M W )F% 138247r 3 (m + M V )F* 

g\Ml i M 2 (2 Cl - c 2 - c 3 ) M 2 (3(J M - Ji 5 ) + 2^ A di 8 ) 



1927r(m + M 7r )m 2 F 2 87i(m + M n )mF 2 6tt (m + M w ) F 2 

2mM 2 (-4e M - 2e 15 - 2ei 9 + e 22 + e 36 ) 



37r(m + M w )F2 

M 2 

- [-36 + ^(133 - 48tt) + 0* (74 + 12tt) 



691 27r 3 (m + M W )F,* 1 

-6m(52ci - c 2 - 32c 3 (l + gfo) + 327rc 4 )] , (A.6) 



mg\ 



24 7 rM^(m + M 7r )F 2 
£4 + 2mc 4 



247r(m + MJF 2 

M w [3g 2 A + 32m 2 (Jx + J 2 ) - 16m 2 g A d 18 ] M v m[Z + ^(21 - 18vr) + g\{2 - 12vr)] 



967rm(m + M W )F 2 34 567r 3 (m + M W )F£ 

g\Ml Mid Ml (Ji + d 2 + 3 J 3 + 2 J 5 + # A d 18 ) 



487r(m + M^)m 2 F 2 I2n (m + M^mF 2 67r(m + M W )F 2 

mM 2 (-2e 17 - 2e 18 - 2e 21 + e 37 ) 



37r(m + MJF 2 
M 2 

-[-18 + g\ (36 + 141tt) + 4(-4 + 42tt) 



69127r 3 (m + M^)F4' 

+8m^i(-127rc 3 + (12tt + ll)c 4 )] , (A.7) 



247rM w (m + A^)F 2 
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9a ~ 4mc 3 
24 7 r(m + M 7r )F2 

M w [3g 2 A + 16mc 2 - 16m 2 (d u - d 15 + g A di 8 )] M^mg\{27>\-K + <d(-56 + 96tt)) 
+ 967rm(m + M„)F% 138247r 2 (m + M V )F% 

g A M 2 d 18 2mM 2 (-4e 14 - 2e 15 - 2e 19 + e 22 + e 36 ) 

48vr(m + M^)m 2 F% 67r(m + M W )F% 37r(m + M n )F% 

M 2 

57— 1 r^[-36 + #1(133 + 24tt) + #i(-10 + 66tt) 

+6m(-52ci + c 2 + (32 - 16tt^)c 3 + 16tt^c 4 )] . (A.8) 
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Figure 1: Fits and predictions for the KA85 phase shifts as a function of the pion laboratory 
momentum q n using strategy one as explained in the text. Fitted in each partial wave are the 
data between 40 and 97 MeV (filled circles). For higher and lower energies, the phases are 
predicted as shown by the solid lines. The other lines refer to the best fits at the various orders 
as explained in the text. 



22 




50 100 150 200 50 100 150 200 



P13 q„[MeV] P1 1 q„ [MeV] 



Figure 2: Fits and predictions for the EM98 phase shifts as a function of q n using strategy 
one as explained in the text. Fitted in each partial wave are the data between 41 and 97 MeV 
(filled circles). For higher and lower energies, the phases are predicted as shown by the solid 
lines. The other lines refer to the best fits at the various orders as explained in the text. 
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Figure 3: Fits and predictions for the SP98 phase shifts as a function of the pion laboratory 
momentum using strategy one as explained in the text. Fitted in each partial wave are the 
data between 60 and 100 MeV (filled circles). For higher and lower energies, the phases are 
predicted as shown by the solid lines. The other lines refer to the best fits at the various orders 
as explained in the text. 
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Figure 4: Fits and predictions for the KA85 phase shifts as a function of the pion laboratory 
momentum q n using strategy two as explained in the text. Fitted in each partial wave are the 
data between 40 and 97 MeV (filled circles). For higher and lower energies, the phases are 
predicted. The various lines refer to the contributions from the various orders as explained in 
the text. 
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Figure 5: Fits and predictions for the EM98 phase shifts as a function of q w using strategy 
two as explained in the text. Fitted in each partial wave are the data between 41 and 97 MeV 
(filled circles). For higher and lower energies, the phases are predicted as shown by the solid 
lines. The various lines refer to the contributions from the various orders as explained in the 
text. 
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Figure 6: Fits and predictions for the SP98 phase shifts as a function of the pion laboratory 
momentum using strategy two as explained in the text. Fitted in each partial wave are the 
data between 60 and 100 MeV (filled circles). For higher and lower energies, the phases are 
predicted as shown by the solid lines. The various lines refer to the contributions from the 
various orders as explained in the text. 
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